Anomalous corrections to the Kelvin equation at complete wetting 
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We consider an Ising model confined in an L x 00 geometry with identical surface fields at the 
boundaries. According to the Kelvin equation the bulk coexistence field scales as 1/L for large 
L; thermodynamics and scaling arguments predict higher order corrections of the type 1/L 2 and 
I/L 5 / 3 at partial and complete wetting respectively. Our numerical results, obtained by density- 
matrix renormalization techniques for systems of widths up to L — 144, are in agreement with a 
1/L 2 correction in the partial wet regime. However at complete wetting we find a large range of 
surface fields and temperatures with a correction to scaling of type 1 / L 4 ' 3 . We show that this term 
is generated by a thin wetting layer whose free energy is dominated by the contacts with the wall. 
For L sufficiently large we expect a crossover to a I/L 5 / 3 correction as predicted by the theory for 
a thick wetting layer. This crossover is indeed found in a solid-on-solid model which provides a 
simplified description of the wetting layer and allows the study of much larger systems than those 
available in a density-matrix renormalization calculation. 
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I. INTRODUCTION 

Wetting phenomena are very common in natureB; be- 
sides the most familiar situation of a liquid- vapor or of a 
binary fluid mixture in contact with a solid wall, wetting 
has also been studied—for instance, at grain boundaries^, 
in polymer mixturesa, alloys 4 !, superconductors^ and 
metallic vaporsO. Usually one considers a semi-infinite 
geometry with a solid planar wall that preferentially ad- 
sorbs one of the phases of a system in thermodynamic 
equilibrium. Below T c , the bulk critical temperature, the 
adsorbed phase forms cither isolated droplets or a thick 
macroscopic layer. The first case, known as partial wet- 
ting, occurs for temperatures below a temperature T w , 
the wetting temperature, while the second case occurs 
for T w < T < T c and it is referred to as complete wet- 
ting. 

In this paper we analyze the effect of wetting on the 
thermodynamics of a two-dimensional Ising model con- 
fined in an L x 00 geometry, with identical surface fields 
at the boundaries. If L — > 00, i.e. in the bulk, phase 
coexistence occurs for temperatures T < T c and for van- 
ishing bulk magnetic field h — 0. It is well known that the 
combined effect of surface fields and confinement shifts 
phase coexistence to a finite value of the bulk magnetic 
field h = h (L) ^ 0, which for large L scales as: 



1 CT COS< 

h o( L ) = - 

L nri}, 



(1) 



where cr , toj, and 9 are the surface tension, the bulk 
spontaneous magnetization and the contact angle, re- 
spectively. Equation (Q) is known in the literature as the 
Kelvin equation_and has a long history that dates back to 
the last centuryQ; it essentially states that to have equi- 
librium between two phases, or phase coexistence, one 



needs a bulk field of the order of 1/L that compensates 
the effect of the surface fields. 

Albano et alu and Parry and Evansi analyzed the next 
order correction term to the Kelvin equation. Both stud- 
ies, using scaling and thermodynamics arguments, con- 
cluded that for temperatures below the wetting temper- 
ature T w , i.e. at partial wetting, the leading correction 
to scaling term is of type 1/L 2 . In the case of complete 
wetting the correction is expected to be non-analytic due 
to a singularity of the s«pface free energy. For the two- 
dimensional Ising modelau the predicted correction term 
is proportional to I/L 5 / 3 . Monte Carlo simulations on 
M x L lattices with M L were also performed!; the 
dominant 1/L scaling of the Kelvin equation w^s well 
verified (this was done also in three dimensions!!!!) , but 
the data were not accurate enough to convincingly test 
the type of corrections to scaling. 

We have studied this, problem using density-matrix 
renormalization groupM (DMRG) techniques and we 
have calculated the bulk coexistence field h (L) for sev- 
eral values of surface fields and temperatures for strips of 
widths up to L = 144. Our numerical data are in good 
agreement with corrections of type 1 /L 2 in the partially 
wet regime, but not compatible with I/L 5 / 3 corrections 
at complete wetting. We find instead that the correction 
to scaling at complete wetting is in good agreement, for 
a wide range of surface fields and temperatures, with a 
term of type 1/L 4 / 3 . An explanation of this discrepancy 
is discussed in detail in the paper. 

The DMRG has been applied mostly to the study 
of pound state properties_pf quantum systems both in 
one!! 2 ! and two dimensions!! 3 ! The method was originally 
developed by WhiteEH for the diagonalization of quantum 
Hamiltonians of very large systems using a truncated ba- 
sis set. It does not have the minus sign problems that 
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plague quantum Monte Carlo simulations of fermionic 
systems at-low temperatures. 

Nishinoli-3 applied White's DMRG to two-dimensional 
classical systems, which are deeply related to quantum 
systems of one dimension lessE-3. In the classical case 
a transfer matrix of a large strip is generated by a se- 
ries of iterations that, starting from a small strip, en- 
large the system and truncate the configurational space 
very efficiently. Although there are obviously no mi- 
nus sign problems in the classical case the advantages 
of the DMRG with respect to the Monte Carlo method 
are mainly the high accuracy and the possibility of deal- 
ing with very large systems. Untill now the method has 
been applied-. only exclusively to two-dimensional classi- 
cal systemal3; although extensions to three dimensions 
are possible, they will require a high computational ef- 
fort and one should not expect the same high accuracy 
that is typically found in two dimensions. 

Recently a series of classical models in a confined ge- 
ometry have been investigated using DMRG techniques; 
the studies-focused on the effect of gravity on phase 
coexistencetH and on the critical point shiftEa for a con- 
fined fluid, on the calculation of critical (density profiles 
and Casimir amplitudes for Potts modelstj and on stud- 
ies of critical adsorptioncd. Another advantage of the 
DMRG for these types of models is that, being based on 
a transfer matrix approach, it naturally deals with ge- 
ometries which are confined along one direction only and 
infinite in the other, i.e. with Lxoo lattices. More details 
on the DMRG can be found in the broad existing liter- 
ature; a short review on the method has been presented 
in Ref. |lj 

The paper is organized as follows: In Sec. || we intro- 
duce the model and present the numerical results focus- 
ing on the corrections to scaling for the Kelvin equation; 
in Sec. Ill we review the existing theory leading to a 
correction to scaling term proportional to 1/L 5 / 3 for the 
two-dimensional Ising model. In Sec. ^ we discuss a 
mechanism that explains the apparent discrepancy of the 
DMRG results with the theory. In Sec. |v| we consider a 
discrete solid-on-solid model and analyze the correction 
to scaling exponent for that case. Sec. terminates our 
paper with some conclusions. 



II. THE DMRG RESULTS 

We consider the Ising model on an L x 00 lattice with 
the following Hamiltonian: 

H = - J } ^SjSj -hi} j s k -hi2_^ s i-h2_^,Si, (2) 

(ij) k I i 

with J > and Si = ±1. The first sum is restricted to 
nearest neighbors, while (I) and (II) indicate the sums 
performed on spins at the boundary only; hi is a sur- 
face field which models the preferential adsorption of 
the boundaries for one of the two phases, and h is a 



bulk magnetic field. In the rest of the paper we take 
J = 1, which fixes the bulk critical temperature to 
T c = l/ln(l + y/2) ~ 2.269. 

The phase diagram for the model for fixed L and hi is 
shown in Fig. 0(a). The thick dashed line indicates the 
bulk phase coexistence line h = 0, T < T c terminating 
at the bulk critical point. For finite L there are no true 
thermodynamic singularities and quantities such as the 
free energy and the magnetization depend smoothly on 
temperature, surface and bulk fields. Fig. |l|(b) shows a 
plot of the total free energy f(h) at fixed T, L, hi: For 
large h, either positive or negative, the system is "frozen" 
with all spins pointing up or down and the free energy as 
a function of h is a straight line with slopes ±1. In three 
or higher dimensions these two branches meet each other 
forming a cusp at phase coexistence or with higher order 
singularities at the finite system critical point. In two di- 
mensions the singularities are replaced by a smooth max- 
imum as indicated in Fig. ^(b). For T < T c we identify 
this maximum with the pseudo-coexistence field ho(L) 
which is expected to scale according to the Kelvin equa- 
tion (Eh. ' 




' T = const. 




FIG. 1. Due to the finiteness of L and to the surface fields 
the pseudo phase-coexistence (thin solid line) is shifted with 
respect to the bulk coexistence line (thick dashed line). The 
temperatures T\ and T2 are in the partially and completely 
wet regimes respectively, (b) The pseudo-coexistence point 
(indicated by the arrow) is identified with the maximum of 
the free energy as function of the bulk field h. 

The free energy maximum at T = T c allows us to 
identify a critical field h c (L), whose scaling properties 
are governed by critical exponents. Finite size scaling 
predicts^ that for large L: 



h c {L) ~ L 



-yn 



(3) 



with j/h = 15/8 for the two-dimensional Ising model. 
Notice that even for T > T c one can find free energy 
maxima; the locus of these points is indicated by a thin 
dashed line in Fig. |l](a) and approaches the h = axis 
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exponentially in L. For our purposes this is an uninter- 
esting region and it has not been investigated. It should 
be stressed that the definition of h c (L) is somewhat ar- 
bitrary. In higher dimensions one can find unambiguosly 
a finite system critical point and it is well known that, 
due to the effect of confinement, both magnetic field and 
temperature are shiftedEa with respect to the bulk criti- 
cal point. Also in the present case it would be possible 
to define a finite system pseudo-critical point identified 
as the point for which, for instance, the specific heat has 
a maximum along the thin solid line of Fig. [l](a); how- 
ever we focus here on the shift along the magnetic field 
direction and the definition given above of h c (L) is ap- 
propriate and definitely less computationally demanding. 

In the rest of this Section we will present the DMRG re- 
sults concerning the scaling analysis of the pseudo-critical 
(h c (L)) and the pseudo-coexistence (ho(L)) fields. We 
start from the first case, which will be just briefly dis- 
cussed. 



A. T = T c 

Figured shows a plot of \n[h c (L)} vs. InL from L = 8 
up to L — 104 and for three different values of the surface 
field hi = -0.1, -0.5 and -0.99. The dotted line has a 
slope equal —15/8, which is the expected value of the 
exponent — j/h for the two-dimensional Ising model and 
is in good agreement with the data only for the highest 
surface fields considered (hi = —0.5, —0.99). 
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of the surface free energy depends on the surface field 
through the scaling variables hxL Al ' 1/ , where Ai is the 
surface gap exponent (recall Ai = 1/2 for the two- 
dimensional Ising model^-. Thus the surface field enters 
in the analysis as followsca: 



(4) 



where K (u>) is a scaling function. For large values of its 
argument the scaling function K(u>) should "saturate" 
to a constant non- vanishing value Koo ^ 0, since one ex- 
pects a scaling behavior as described in Eq. @j for suf- 
ficiently large L. In absence of surface fields there is no 
shift along the magnetic direction, i.e. h c (L) — 0, which 
implies K(0) = 0. For small ui one expects K(u) ~ ui 
since h c (L) should scale linearly with the surface field. 
This leads to the following prediction: 



h c (L) ~ L 



-(i/h-Ai/i/) 



(5) 



for hi small and L not too large. This scaling behavior 
is clearly seen in the numerical DMRG data for h c (L) 
shown in Fig. || for the lowest surface field analyzed, i.e. 
h\ = —0.1, where the dashed line has a slope equal —11/8 
in agreement with the prediction of Eq. (^|). From the 
data for hi = —0.5 one finds an estimate of the crossover 
size at L « 40. The shape of the scaling function K(ui) 
found from the DMRG data is shown in the inset (a) of 
Fig. U and it is indeed linear for small u> and it is found 
to saturate to a constant ss 10 for u) > loq ~ 5. 

To verify that our numerical data fulfill a scaling of the 
type presented in Eq. we have plotted in the inset 
(b) of Fig. U again ln[/i c (L)] vs. InL for constant values 
of the scaling variable u = hiL Al ^ . The data follow 
nice straight lines with slope —15/8 even in the region of 
small u>, i.e. for uj <C wq. 

Finally our estimate of the magnetic exponent j/h 
from the data of h c (L) at hi = —0.99 and —0.5 yield 
yu = 1.874999(1) and y H = 1.87506(1), respectively, in 
an excellent agreement with the exact value 15/8. Such a 
high numerical accuracy for critical exponents is typical 
of the DMRG method used in combination with powerful 
extrapolation techniquesE^ . 



B. T < T c 



InL 



FIG. 2. Crossover between an exponent yH — Aij/t = 11/8 
(the dashed lines) and an exponent yu = 15/8 (the dotted 
lines) for the scaling of the critical field h c (L). Insets: (a) 
Plot of the scaling function K(u>) vs. ui calculated from the 
numerical DMRG data, (b) Plots ln[/i c (L)] vs. InL for con- 
stant values of the scaling variable lj = h\L AlVT . 

The behavior of h c (L) for smaller surface fields can 
be understood from finite-size scaling: The singular part 



For temperatures below the bulk critical temperature 
T c ~ 2.269 one expects a scaling of the type ([!]). In 
practice however, since only a limited size L is available 
for numerical computation, it is preferable to consider 
only temperatures not too close to T c , where a crossover 
towards a scaling of the type of Eq. (||) is expected. At 
low temperatures one also has the advantage that DMRG 
calculations become very rapid. In the present study we 
have calculated a series of values of ho(L) up to L = 144 
and for temperatures ranging from T ~ 0.97 T c down to 
T w 0.4 T c with different values of the surface field hi- 
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We introduce the logarithmic derivatives: 
ln[/i Q (i)]-ln[/io(i + 8)] 



XL 

for L = 8, 16, 24, 
pansion: 



(6) 



ln£-ln(i + 8) 
. , 136. Assuming the following ex- 



w = 4 



£7 



(7) 



where we expect a = 1 and A given by the Kelvin equa- 
tion (Q) and 7 > a a correction-to-scaling exponent, one 
has to lowest orders in 1/L: 



XL 



a 1 



1 



B 



A L-y- 



(8) 



Figure || shows a plot of xl vs. 1 / L for different values 
of temperatures and surface fields. The set (a) corre- 
sponds to the case of partial wetting, where the scaling 
behavior oi xl should be linear in 1/L since one expects 
7 = 2 and a = 1 in Eq. (g). The fact that these data 
follow straight lines confirms the behavior predicted by 
the theory. 

A clearly different behavior is observed for complete 
wetting, as shown in the sets (b) and (c) and in the inset 
of Fig. g. In this case Xl has a non-monotonic behavior 
as a function of 1/L; the numerical data do not follow 
straight lines and seem to approach the value a = 1 as 
1/L P , with p < 1. Notice the difference between the 
data at low temperatures (sets (b) and (c)) and at tem- 
peratures not too far from the bulk critical temperature 
(inset). In the former case the value of the logarithmic 



derivative xl for the largest size analyzed is just a few 
percent off the expected asymptotic value Xl — * ot = 1. 
In the latter case Xl is much further away from its ex- 
pected asymptotic value. This is the effect of a crossover 
towards a scaling of the type (||) governed by the mag- 
netic exponent j/h = 15/8, which is shown by a cross in 
the inset of Fig. 0. 
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FIG. 3. The plots of xl vs. 1/L for different values of tem- 
perature and surface field where (a) is in the partially wet 
regime, while (b), (c) and the curves in the inset are in the 
completely wet regime. The value of the surface field for the 
inset is hi — —1.0; notice the trend to a crossover towards 
the magnetic exponent j/h = 15/8 indicated by an arrow in 
the inset. 



TABLE I. The values of the scaling exponent a and the coefficient A obtained by the BST extrapolationl-3 of the logarithmic 
derivatives xl in the limit L —* oo, for different values of T and hi. The third column shows the values of the wetting 
temperature T w (hi), which is known exacttytJ in the two-dimensional Ising model. In the last column we show the ratio of 
surface tension cro and bulk magnetization ms calculated from the exact solution of the two-dimensional Ising model, which is 
expected to be identical to the coefficient A (see Eq. [j]) in the complete wet regime. The convergence to the expected value 
a = 1 is always very good, but it is clearly better for partial wetting, as compared to complete wetting. 
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hi 


T w (hi) 
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<Jo/m b 


1.0 


-0.5 


1.958... 


1.000000(1) 


0.99077(1) 


1.72891 


1.0 


-0.927... 


1.0 


1.0003(1) 


1.729(1) 


1.72891 


1.0 


-1.0 





1.0003(1) 


1.729(1) 


1.72891 


2.1 


-0.5 


1.958... 


1.004(2) 


0.335(2) 


0.33508 



Table |l| shows the values of a and A obtained from an 
accurate numerical extrapolation of the DMRG data for 
different values of the temperature and surface field, both 
at partial (T < T w ) and at complete (T w < T < T c ) wet- 
ting. For the extrapolation we used the BST algorithmcJ. 
Notice the excellent agreement with a = 1 as predicted 
by the Kelvin equation (Q) . Extrapolations are better at 
low temperatures and at partial wetting, in agreement 
with the fact that higher order corrections in 1/L vanish 



more rapidly in the limit L — > oo with respect to the 
complete wetting case. 

The last column of Table | shows the ratio between the 
surface tension (cr ) and bulk spontaneous magnetization 
(nib), which are known exactly in the two-dimensional 
Ising model. According to the Kelvin equation (]l|) this 
ratio should be equal to the coefficient A at complete 
wetting where cos 9 — 1. The numerical results agree 
very well with this; notice also that at complete wetting 



4 



the coefficient A does not depend on tfie surface field 
anymore. 

Having found a — 1 with high accuracy we can now 
consider higher order corrections; setting a = 1 in Eq. 
(H) one obtains: 



In (x l — 1) = In 



( 7 -l) lnZ + 



(9) 



Figure U shows a plot of In (xl — 1) vs. InL for different 
temperatures and surface fields. The behavior at par- 
tial (a) and complete wetting (b,c and d) can be clearly 
distinguished. The straight dashed line in the figure has 
slope —1 and fits well the asymptotic regime at partial 
wetting, in good agreement with the expected correction- 
to-scaling exponent 7 = 2. Other calculations at partial 
wetting are in good agreement with this values as well. 




FIG. 4. Plots of ln(x_L — 1) vs. InL for different values of 
surface field and temperature. The different sets correspond 
to (a) T = 1.0, hi = -0.5 (b) T = 1.0, hi = h w = -0.9271 (c) 
T = 1.0, hi = -1.0 (d) T = 2.1, hi = -0.5 and (e) T = 2.2, 
hi = —1.0. The dashed line has slope —1, whereas the dotted 
lines have slope —1/3. 

The asymptotic behavior for the sets (b), (c) and (d) 
of Fig. H is well-fitted by straight lines in the log-log plot 
with slope —1/3 (dotted lines). This is not in agree- 
ment with the current theory which predicts a correction- 
to-scaling exponent 7 = 5/3 at complete wetting and 
one would therefore expect an asymptotic slope equal to 
1 — 7 = —2/3, i.e. twice as large than that of the dot- 
ted lines of Fig. ^. For the set (e) one observes a slight 
deviation of the data from —1/3 towards higher slopes 
at large L, but not compatible with —2/3. This suggests 
that the large value of L considered here (L = 144) is 
still not sufficient to be in the asymptotic regime. The 
striking characteristic of the data at complete wetting is 
that for a large range of temperatures and surface fields 
(see sets (b), (c) and (d)) and up to L = 144 the cor- 
rection to scaling seems to be given by a term of type 
I/L 4 / 3 . This behavior will be discussed in more detail in 



the next sections. 

To complete the presentation of the DMRG results we 
show in Fig. || the magnetization profiles for the two 
coexisting states. The magnetizations are plotted as a 
function of the scaled variable l/L, with I = 1, 2, . . . , L 
labeling the distance from one of the walls. Temperature 
and surface fields are fixed to T = 2.1 and hi = —0.5 
respectively; this corresponds to a regime of complete 
wetting (T w (hi = —0.5) = 1.958). We should also stress 
that the profiles refer to bulk fields slightly lower and 
higher than the pseudo-coexistence field ho(L). For bulk 
field exactly equal to ho(L) the magnetizations of the two 
coexisting phases are averaged and it is not possible to 
distinguish between them. The positive bulk field favors 
a bulk phase with positive magnetization but the nega- 
tive surface fields favor the adsorption of negative spins 
at the boundaries. Since T > T w the negative spins form 
a layer that wets the walls (as can be seen clearly in 
the three uppermost profiles of Fig. [5]) and thickens as 
L — * 00, or h — ► 0. As will be derived in the next section 
the wetting layer grows fofi large L as L 1 / 3 , thus it re- 
mains macroscopically thinEJ in the thermodynamic limit 
L — > 00. This can be seen in the figure as well, since 
in the scaled plots of Fig. || the wetting layer flattens 
towards the wall as L increases. 
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FIG. 5. The magnetization profiles of the two coexisting 
phases at T = 2.1 and hi = —0.5, approaching complete 
wetting, for different strip widths (circles: L = 150, squares: 
L = 300 and triangles: L = 500). 



III. THERMODYNAMIC THEORY 

We briefly review here the theoretical arguments lead- 
ing to a correction-to-scaling exponent 7 = 5/3 for the 
two-dimensional Ising jpodel, following the derivation 
given by Albano et aln. A similar derivation, but ex- 
tended to higher dimensionalities and also to long range 
interaction forces has been presented in Ref. ^[ 
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For T < T c the system can in be in one of the two 
phases with the majority of spins pointing either up or 
down with magnetization profiles similar to those shown 
in Fig. [|. The total free energy can be divided into bulk 
and surface contributions as follows: 



f ± (h,h u T) = f±(h,T) + -f™±(h,h 1 ,T) 



(10) 



where the superscripts + or — indicate the two phases 
and the surface term f w± is the excess free energy for 
the wall-spin up (or down) phase. There is no interaction 
between the surfaces since, as pointed out in the previous 
section, even for T > T w the surface layers remain very 
far from each other in the limit L — > oo. 

Following Ref. || we expand the bulk term to lowest 
orders in h: 



f b ± (h,T) = f b (T) T m b h + ^h 2 



(11) 



with mb > the bulk spontaneous magnetization and x 
the susceptibility. For the surface term we write: 



f^(h,h!,T) 



+ mf(h, hi, T) h 



(12) 



where a w± — ff ± (0, h\, T) and mf is the surface excess 
magnetization. The condition for phase coexistence is: 



f + (h,h 1 ,T) = f-(h,h 1 ,T) 



(13) 



From this relation and from Eqs. (|il| ) and (|l2|) one finds 
to lowest orders in 1/L the following expression for the 
coexistence bulk field: 



ho(L) = \ 



w + 



(14) 



m b - ~ (m^ - m s ) 

Using Young's relation^ a w + — a w - — <7o cos#, one finds 
to leading order in 1 /L the Kelvin equation ([!]) . 

At partial wetting the surface magnetizations re- 
main finite and higher order corrections to the Kelvin 
equation are expected to be of type in agreement 

with our numerical results. More interesting is the case 
of complete wetting where the thickness of the two wet- 
ting layers in the two-dimensional Ising model diverge 
as I ~ ft -1 / 3 for h — > 0. From h ~ 1/L one derives 
I ~ L 1 / 3 , which is the result anticipated at the end of 
the previous section. Since the surface magnetization di- 
verges with the same power of h as the thickunss I, i.e. 
m s ~ ft -1 / 3 ~ L 1 / 3 , from Eq. ( |l4| ) one obtainstrQ: 



ho(L) 



A 
L 



B 

IsT 3 



(15) 



I is large enough for the scaling relation m s ~ I ~ /i -1 / 3 
to be valid. For thin layers one expects a deviation from 
this behavior and consequently from a correction to scal- 
ing of the type 1/L 5 / 3 . Figures ||(a) and ||(b) show two 
possible configurations: A very thin wetting layer with 
many contacts with the wall (a) and a thicker one with 
rare contacts and with an interface well-separated from 
the wall (b) . Although the DMRG calculations were per- 
formed up to large L it seems quite clear from the mag- 
netization profiles of Fig. || that the typical configuration 
of the wetting layer is closer to that of Fig. ^|(a) than to 
that shown in Fig. ||(b). In the case of a thick wetting 
layer, as in Fig. ^(b), one would expect magnetization 
profiles that start off with a plateau close to the wall and 
decay to the bulk magnetization quite far from it. No 
signs of a plateau are to be seen in the profiles of Fig. ||: 
The decay to the bulk magnetization is very rapid and 
starts already from the boundary spins. 



(a) 



(b) 




IV. THICKNESS OF THE WETTING LAYER 



FIG. 6. (a,b) The two regimes for approach to complete 
wetting. In the case (a) the contacts with the wall produce a 
surface free energy proportional to h 1 / 3 which results in the 
correction to the Kelvin equation of the type 1/L 4//3 . In the 
case (b) when the interface is sufficiently far from the wall 
this free energy is proportional to h 2 ^ 3 and the correction 
to the Kelvin equation, according to existing theoriesQ'B, are 
expected to be of type 1/L 5 ^ 3 . (c) The dependence of the sin- 
gular part of the surface free energy on the field; the vertical 
dashed line separate the two different scaling regions. 



To explain the corrections to the Kelvin equation found 
in the DMRG calculations we have to understand the 
scaling properties of the singular part of the free energy 
of a thin wetting layer, such as that shown in Fig. [|(a). 
One expects a free energy proportional to 1/h, i.e. pro- 
portional to the probability of contacts with the wall (l\ 
is the average distance between two contact points) . The 
scaling behavior of l\ is well known in the literature (see, 
for instance, Refs. |26| , |2~7| ) and it is related to the mean 
first length of a return to the "origin" , i.e. the wall, for 
a random walk which is biased towards the wall due to 
the presence of the bulk field. It has been shown that: 



The argument leading to the expansion ( |l5| ) is based 
on the assumption that the thickness of the wetting layer 



h 



-1/3 



(16) 



G 



Assuming thus that the surface free energy U) has a sin- 
gular term proportional to ft, 1 / 3 and repeating the analy- 
sis of the previous section one finds a correction-to-scaling 
term of the type 1/L 4 / 3 . The case of Fig. ||(b) has been 
treated in the previous section; when I, the thickness of 
the wetting film, is large enough the singular part of the 
surface free energy is dominated by the magnetic energy 
of this layer and thus proportional to m s h ~ Ih ~ /i 2 / 3 . 

This simple scaling theory predicts a crossover from 
thin to thick wetting layer which results in the differ- 
ent singular behavior of the surface free energy shown in 
Fig. ^(c). As a consequence one expects a crossover in 
the correction-to-scaling term from 1/L 4 / 3 to 1/L 5 / 3 . 



V. THE SOLID-ON-SOLID MODEL 

In order to obtain a thicker wetting layer and observe 
the predicted crossover one could increase the tempera- 
ture in the DMRG calculations, since entropy favors an 
interface far from the wall where more configurations are 
available. This is indeed consistent with the behavior of 
the data of Fig. |]; at the largest temperature considered 
(the set (e) corresponds to T « 0.97 T c ) a deviation from 
a correction to scaling of the type 1 / L 4 / 3 towards higher 
powers of 1/L can be clearly observed. This deviation is 
not yet compatible with the expected 1/L 5 / 3 . However 
one cannot increase the temperature too much since close 
to T c there is a crossover to a scaling of the type given in 
Eq. (H), which makes the analysis very difficult. 

A second possibility is to enlarge the strip width. Al- 
though DMRG calculations for L > 144 are feasible (we 
have calculated magnetization profiles up to L = 500 as 
shown in Fig. ||) it turns out that it is difficult to keep 
good accuracy for larger L since Hq (L) decreases and one 
needs higher computational effort to maintain a small 
relative error on the data. 

In this section we analyze the correction to scaling 
to the_Kelvin equation (0) in the solid-on-solid (SOS) 
modelEj which provides a simplified description of the 
wetting layer. The model does not have a bulk critical 
point (thus no crossover to a scaling of type @ is ex- 
pected) and it allows the study of much larger systems 
than those available in a DMRG calculation. 

In the SOS model a sharp interface is assumed to sep- 
arate the two phases with all spins pointing either up or 
down; fluctuations within each phase are neglected. No 
overhangs are allowed, so the position of the interface is 
described by a single valued discrete function l(k), where 
k labels the coordinate position along the wall and I is 
the distance of the interface from the wall. 

Neglecting the interaction between the wetting layers 
at the two walls we take a single wetting layer in a strip of 
size L/2, so the discrete coordinate takes the values 1 = 0, 
1, . . . , L/2. We consider the one-dimensional transfer 
matrix, whose elements are given by Tw = exp(— (3£u>), 
where (3 is the inverse temperature and the energy is 



given b 



E w =2J\l-l'\ + (J + h 1 ) (0,+, 
+ h(l + l' - 1 



with: 



if 1 = 

1 if I > 



(17) 



(18) 



The configuration l(k) = corresponds to all spins up 
and gives the ground state energy: 



eo 



-hi 



hL 
~2~' 



The free energy per row is equal to: 



(19) 



(20) 



with A max the largest eigenvalue of the transfer matrix 
T. For the other phase where no wetting layer is present 
we assume that all the spins are pointing down, which 
gives the free energy: 



f sos (h,hi,T,L) = hi 



hL 



(21) 



The value of the bulk field which gives phase coexistence, 
ho(L), can be found by equating the two free energies in 
Eqs. ( po|) and (|2l]), which we have done numerically for 
strips of widths up to L = 2000. The advantage of the 
SOS approach is that it allows one to deal with very large 
systems since the size of the transfer matrix T grows only 
linearly as function of the strip width L. 

Figure [7] shows a plot of \u(xl — 1) vs. InL, where 
xl is given in Eq. (^), for various values of temperature 
and surface field. The sets (a) and (b) refer to partial 
wetting and their asymptotic behavior is very well fitted 
by straight lines with slope —1, in agreement with a cor- 
rection to scaling of type 1/L 2 for this case. The other 
sets refer to T > T w , i.e. to complete wetting. The set 
(c) corresponds to T = 1 and hi = —0.95. For a wide 
range of system sizes (56 < L < 250) the SOS data fol- 
low closely a straight line with slope —1/3, as also found 
in the DMRG calculation, but for L large enough the 
data bends towards higher slopes, demonstrating that 
the correction to scaling found in DMRG calculations is 
not the true asymptotic behavior. The crossover region 
appears to be quite large since still for L = 2000 we do 
not find the expected slope —2/3. At this low tempera- 
ture T = 1 » 0.4 T c the SOS data are also in quantitative 
good agreement with the DMRG results. At higher tem- 
peratures this is not the case anymore. The set (d) refers 
to T = 2.1, hi = —1.0. In this case there is a good agree- 
ment with the slope —2/3 of the dashed-dotted line, as 
expected on the basis of the scaling theory of Albano et 
al.u and of Parry and Evansa. It is also important to 
point out that for the full Ising calculation at the same 
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temperature and surface field we still find a slope —1/3, 
which suggests that the Ising interface is closer to the 
wall than the SOS interface for the same values of sur- 
face field and temperature. 

We have also calculated (I), the average thickness of 
the wetting layer, from the SOS Hamiltonian. For the 
sets (c) and (d) and L = 200 we found (I) = 1.3 and 
(I) = 11.5 respectively (these values are expressed in 
terms of lattice units). Notice that for a system of size 
L = 200 one observes corrections to scaling of the type 
1/L 4 / 3 and 1/L 5 / 3 in the sets (c) and (d) respectively and 
the large differencee in the values of (I) is consistent with 
a crossover from a thin to a thick wetting layer discussed 
in the previous section. 
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FIG. 7. Plot of \n(xL — 1) vs. InL in the SOS approxima- 
tion for (a) T = 1.0, hi = -0.3, (b) T = 1.0, hi = -0.7, (c) 
T = 1.0, hi = -0.95 and (d) T = 2.1, hi = -1.0. The set (a) 
and (b) correspond to the partial wet regime, while (c) and (d) 
to complete wetting. The slopes of the straight lines is equal 
to —1 (dashed), —1/3 (dotted) and —2/3 (dashed-dotted). 



VI. CONCLUSIONS 

In this paper we have analyzed the influence of wet- 
ting on the thermodynamics of an Ising model in a con- 
fined geometry with equal surface fields at the bound- 
aries. Firstly, we verified with high numerical accuracy 
that the coexistence bulk field scales as 1/L (L is the 
distance between the boundaries) in agreement with the 
predictions of the Kelvin equation. Then our study fo- 
cused on the higher order corrections to this equation 
which are expected to be analytic (non-analytic) at par- 
tial (complete) wetting. 

Using density-matrix renormalization techniques we 
found that for a large range of surface fields and temper- 
ature higher order corrections are not compatible with 
1/L 5 / 3 , which is the behavior predicted by the existing 
theory, but they are of type 1/L 4 / 3 . We have shown that 



this apparent disagreement is due to the fact that even 
for the large sizes considered (L = 144) the wetting layer 
has a limited thickness, so that the singular part of the 
surface free energy which determines the correction-to- 
scaling behavior is dominated by the contacts with the 
walls. In this case a simple random walk argument in- 
deed predicts a correction term of the type 1 /L 4 / 3 for a 
thin wetting layer. We have analyzed the problem also 
using a solid-on-solid model, which provides a simplified 
description of the wetting layer, but allows the study of 
much larger systems (strips of widths up to L = 2000 
were considered). The results of this analysis are consis- 
tent with the 1/L 4 / 3 correction for small thicknesses of 
the wetting layer, but also show a crossover to a correc- 
tion of type 1/L 5 / 3 for sufficiently thick layers. 
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